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Two symmetries of the local nonlinear Sf gyrokinetic system of equations in tokamaks in the 
high flow regime are presented. The turbulent transport of toroidal angular momentum changes 
sign under an up-down reflection of the tokamak and a sign change of both the rotation and the 
rotation shear. Thus, the turbulent transport of toroidal angular momentum must vanish for up- 
down symmetric tokamaks in the absence of both rotation and rotation shear. This has important 
implications for the modeling of spontaneous rotation. 
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I. INTRODUCTION 



Understanding how momentum is distributed within 
a tokamak is crucial to explaining experimentally ob- 
served regions of reduced turbulent transport like In- 
ternal Transport Barriers (ITBs) [1] and edge pedestals 
[2]. In these regions of reduced transport, turbulence is 
suppressed by strong velocity gradients, which are deter- 
mined by the transport of toroidal angular momentum. 
Nonlinear Sf flux tube simulations |3j-{6j are extensively 
used to study turbulent transport of toroidal angular mo- 
mentum in tokamaks [THTl]. The time averaged turbu- 
lent transport of momentum is found to be zero for up- 
down symmetric tokamaks in the absence of toroidal ve- 
locity and toroidal velocity gradient. This behavior was 
addressed in Ref. p^, where the transport of parallel 
momentum due to electrostatic fluctuations was studied 
with a quasilinear model. Ref. [15] showed that the lin- 
earized gyrokinetic equation has a symmetry that results 
in the transport of parallel momentum vanishing for up- 
down symmetric tokamaks in the absence of velocity and 
velocity shear. Moreover, this analysis was successfully 
employed to identify the up-down asymmetry as a drive 
for momentum transport [5]. 

However, there are six major effects that were not 
treated in the analysis of Ref. |15j . First, the symmetry 
given in Ref. [15] is only valid for the linearized equa- 
tion. Second, the formulation was collisionless. Third, 
the effect of the E x B shear on transport of momentum 
[IOllTl] [T6] was not considered. Fourth, the Mach number 
was assumed to be small and the effect of the centrifugal 
force [12] was neglected. Fifth, in the evaluation of the 
radial transport of toroidal angular momentum, the con- 
tribution of the component of the velocity perpendicular 
to the magnetic field was neglected. Sixth, electromag- 
netic effects were not taken into account. In this article, 
we present the complete symmetry argument. We con- 
sider the high flow regime, for which the mean velocity is 
comparable to the thermal speed of the ions, and we use 



the local approximation, for which the turbulent trans- 
port only depends on local quantities. In this limit, the 
sign of the radial transport of toroidal angular momen- 
tum changes under a sign change of the average velocity 
and its gradient and an up-down reflection. Thus, for 
an up-down symmetric tokamak without any velocity or 
velocity shear, the momentum transport is zero. 

The symmetry presented here is broken if higher or- 
der terms in the gyrokinetic expansion in = Pi/a <^ 1 
are kept, giving transport of momentum that is small 
in e,; compared to the high flow results. Here pi is the 
ion gyroradius and a is the minor radius of the tokamak. 
The effect of these higher order terms on transport of 
momentum is discussed in Ref. [T7], where a higher or- 
der equation valid in the limit of small poloidal magnetic 
field is proposed. The complete higher order gyrokinetic 
Fokker-Planck equation valid for any magnetic geome- 
try was self-consistently calculated for the first time in 
Ref. [18] for the electrostatic approximation. 

The remainder of this article is organized as follows. 
We first present the local 5f gyrokinetic model in the high 
fiow regime in section|Tlj We then show in section [Hi] that 
the gyrokinetic equations have two fundamental symme- 
tries that involve an up-down reflection. Finally, the con- 
sequences of these symmetries for transport of momen- 
tum are discussed in section IIVI The most cumbersome 
parts of the calculation are given in Appendices [^][Gl 



II. LOCAL Sf GYROKINETIC FORMULATION 

In this section we present the local Sf gyrokinetic sys- 
tem of equations in the high flow ordering, for which the 
average velocity of the ions Vj is of the order of the ion 
thermal speed vu- 

First, in subsection II A we review the magnetohydro- 
dynamic (MHD) equilibrium for a tokamak and we dis- 
cuss the background density, temperature and velocity. 



We then present in subsection |IIB| the system of equa- 
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tions for the fluctuating pieces of the distribution func- 
tion and the electromagnetic field. These equations are 
simplified in the local limit in subsection |II C| We finish 
by giving the radial flux of toroidal angular momentum 
as a function of the fluctuating distribution function and 
electromagnetic fields in subsection |II D[ 



A. Equilibrium 

The background axisymmetric magnetic field of a toka- 
mak is given by 



B = /VC + VC X VV', 



(1) 



where "ip is the poloidal magnetic flux, C is the toroidal 
angle, VC — C/R, with ^ the unit vector in the toroidal 
direction and R the major radius, and / = RBi^, with 
the toroidal magnetic field. We use ( and a poloidal 
angle 9 as our coordinates. The Jacobian of the trans- 
formation from cartesian coordinates x into ip, and C, 
is 



J = KW X V6I) • VCr^ = |B • V6i| 



(2) 



where we take the absolute value to avoid having a sign 
in the Jacobian. 

In the high flow ordering, the electric field is electro- 
static and purely radial to lowest order, and the poten- 
tial is given by e$_i('(/', t)/Te , with Tg the electron 
temperature and e the magnitude of the electron charge. 
To lowest order, all species rotate at the same speed 
within a flux surface, i.e., = = Rn^{ip,t)(!^ ^ 
[THHS]) where the subscript s indicates the species, and 
= — c(9$_i/9'0) is the rotation frequency. To 
next order, there is an axisymmetric component of the 
potential, e$o('0j ^)/T'e ~ 1 that ensures quasineutral- 
ity within a flux surface. 

It is convenient to solve the equations in the frame ro- 
tating with r2^. Thus, from now on we use the peculiar 
velocity w = v — where v is the velocity in the lab- 
oratory frame. In the frame rotating with $7^, the lowest 
order distribution functions are stationary Maxwellians 
Jms = ns{ms/2'KTs)^^'^ exp{—msw'^ /2Ts), where nis is 
the species mass, and rig and Tg are the species den- 
sity and temperature, respectively. The temperatures 
Ts{tp,t) are flux functions. Due to the centrifugal force 
and the electrostatic potential, the densities are not con- 
stant in a flux surface. Instead, they are given by 



n^itp, 9, t) = t) exp - 







msQ-lR^ 

2r, 



(3) 



with 7]s{ip,t) a flux function and Z^e the species charge. 
Imposing quasineutrality. 



(4) 



gives the piece of the electrostatic potential ^o{ip,9,t). 

With these lowest order equilibrium distribution func- 
tions, I{'ip,t) = RBq is a flux function, and the poloidal 
component of the magnetic field is determined by the 
Grad-Shafranov equation 



dtp 



AttR' 



dp 



where the pressure is 



pi^ip, 9,t) =2_^ TjsTs exp - 



Zse'^ 
T, 



rrisVilR'^ 
2Ts 



and 



dp 



dp 



dR 
d^' 



(5) 



(6) 



(7) 



B. Equations for the turbulent fluctuations 

The 6f formulation in the high flow regime is discussed 
in P^25| . It determines the turbulent pieces of the dis- 
tribution function, ^/s(x, w,t) ^ eifMs, the electrostatic 
potential, 0(x,t) ~ CiTe/e, the parallel component of 
the vector potential, ^||(x,t) ^ £iPiB, and the parallel 
component of the magnetic field, i3||(x,t) ^ e^S. The 
perturbed electric and magnetic fields in the frame ro- 
tating with velocity fi^ are to lowest order (5E = — V_l0 
and SB = V_Lyl|| x b + -B||b, where b = B/B is the unit 
vector in the direction of the background magnetic field, 
and Vj^ = V — bb- V is the gradient perpendicular to the 
background magnetic field. We assume that the perpen- 
dicular gradient of the fluctuations is large, on the order 
of the inverse of the ion gyroradius, Vj^ ^/Pii whereas 
the parallel gradient is comparable to the inverse of the 
characteristic length of the device, b- V ~ 1/a. The over- 
line ~ over the perturbed quantities will be useful later 
because the quantities without this ovcrline will be their 
Fourier transforms. 

To eliminate the fast gyrofrequency time scale, it is 
necessary to change to gyrokinetic phase space coordi- 
nates. We use the guiding center position X = x -|- 
il^^w x b, the velocity parallel to the magnetic field 
It'll = w • b, the magnetic moment fi ~ w\/2B and the 

gyrophase ip = arctan[(w x b) • \7'il;/{w ■ Vi/')], where 
ils = ZseB/rUsC is the species gyrofrequency and w_l = 
w — W| I b is the velocity perpendicular to the background 
magnetic field. With these new coordinates, it is easy 
to average out the gyrofrequency time scale. Extracting 
the Maxwell-Boltzmann response from the distribution 
function, i.e., 5fs = - {Zse(j)/Ts)fMs, hsO^,w\\, iJ,,t) 
can be proven to be gyrophase independent, and the 6f 
kinetic equation becomes 



dhs 
dt 



Rfl(C ■ Vx/is + W||b • Vx/is + ^d,s ■ Vx/is 
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ZsefMs ( digU 



T, V B 



Vx$o - R^jyyiR 



(8) 



it is convenient to write the collision operator using the 
Rosenbluth potentials [26]. The particular expressions 
that we use are given in Appendix \K\ 

The fluctuations of the electrostatic potential 
are obtained using the quasineutrality equation, 
J2s Zs J d'^w 6fs — 0, where the integral in velocity space 
is taken at fixed x. Employing (Pw — B dw^^ d^dcp, we 
find that quasineutrality becomes 

^ ZgB I dw\\ d^dtp hs{x — ps) — ^ Zgens(j> ^ 



where the Maxwellian is written in the new variables as 
Ims = n,(m,/2^T,)3/2exp[-m,(u;2/2 + /iS)/r,]. The 
drift 



B 



Vx(x)s X b 



(9) 



accounts for the parallel motion along the perturbed 
magnetic field lines and for the turbulent E x B and WB 
drifts. The generalized potential {x)s is given by 

{X)s = (0(X + ps,t))s~ (X + p„ t))s 



(14) 

where we have stressed that hg depends on velocity 
not only through w\\ and /x, but also through X = 
X — ps{fJ-,ip). The parallel vector potential A|| and 
the parallel magnetic field -B|| are both determined 
by the current equation V x (VA|| x b + i?||b) = 
(47r/c)^gZse J d^w SfgW. The parallel component of 
this equation gives 

^2 1 AttZscB /■ , - 
- V_lA|| = 2^ / dwii dfj.d(phs(x - Ps)w\\, 



{wj_ ■ Aj^(X + ps,t))s, (10) and the perpendicular component is 



(15) 



where A_l is related to -B|| by = b • (V x Aj^) and 
(. . .)s is the gyroaverage holding X, w\\, fi and t fixed. 
The dependence of (j), A^^ and Aj_ on the gyrophase is 
through x = X + Ps that contains the gyroradius Ps — 
— r^^^w X b, also given by 



V2fiB 



[— sin(/3 Vt/i + cos(^ (b x Vip)]. (11) 



The perpendicular drifts v^.s are 



(12) 



where v^o = — (c/i?)Vx$o x b is the E x B drift due 
to the background axisymmetric potential $oj 'Vai,s — 
(/i/17s)bx Vx-B + (w|j/17s)bx (b- Vxb) are the VB and 
curvature drifts, Vco.s = {2w\\ili^/ns)h x [(Vx^ x C) ^ b] 

is the Coriolis drift, and \'cf,s — — (i?ri^/f^s)b x Vx^ is 
the centrifugal force drift. The parallel acceleration is 



-^ib - VxB 



b- Vx'I'o + i^f^cb- Vx-R- (13) 



Finally, {cf))s (C^s' [^s(x - Ps, W|| , i), /ms'])s + 
{Css' [fMs , /is' (x - p^/ , W| I , ^, t)]) s is the gyroaveraged lin- 
earized collision operator between species s and s' . Here 
Css'[fsT fs'] is the bilinear full collision operator, and we 
have emphasized that /is (X, wi | , n, t) and /is' (X, W| | , /i, t) 
depend on velocity space not only through the variables 
W|| and ^1 but also through X = x — ps{^, ^p). The lin- 
earized collision operator is linear in both hg and hgi . To 
show the symmetries of the local Sf gyrokinetic model. 



Vi3|| X b = 2, / dw\ \ dfi dip /is(x — ps)wj^. 

(16) 

We have neglected the parallel gradient b-V^l/aofv4|| 
and the gradient V ^ 1/a of the background magnetic 
field B because they are much smaller than the perpen- 
dicular gradients V_l ~ 1/pi of the fluctuations. 



C. Local approximation 

Since the characteristic perpendicular length of the 
turbulent structures is much shorter than the parallel 
length, it is convenient to describe the turbulent pieces 
/is(V'(X), a(X, t),eiX),w\\, /), (^(V-lx), a(x, t), 0(x), t), 
All (V-Cx), a(x, /), 6l(x), t) and ^|| (V'(x), a(x, t), 6l(x), t) by 
two coordinates perpendicular to the magnetic field line, 
■0 and a, and a coordinate to locate the position along 
the magnetic field that in this case is the poloidal angle 6. 
We stress that the distribution function depends on the 
guiding center position X and that the electromagnetic 
fields depend on the position x. Here Q!(x,i) is the co- 
ordinate in the direction perpendicular to the magnetic 
field within the flux surface and rotating with toroidal 
angular velocity f]^, defined such that B — Va x Wtp 
and da/dt + fl^ida/dC) = 0, i.e.. 



C-g(V)^(V',^)-f^c(^K 



where 



m 



/ d9' (R^B-We)-^ 
Jo 



(17) 



(18) 
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and g(V') = \{I/2Tr)J^'' d0{R^B ■ V9)-'^\ is the safety 
factor. The advantage of using the time dependent vari- 
able a is that the E x B shear [11] is included in the 
formulation in a natural way, as we will see. 

In the local approximation, the characteristic perpen- 
dicular size of the turbulent structures is assumed to be 
small compared to the background length scale. The 
local Sf equations only depend on tp and a through 
the unknowns hs, (f>, ^|| and -B||. It is then more 
convenient to Fourier analyze the distribution func- 
tion and the electromagnetic fields. Defining g = 
(A^Aa)"^/^^ dip J^^ da ge^p{-ik^ip-ikaa), with i = 

V^, 9 = 9 = and A^ 

and Aa the size of the domain of interest, and using Ap- 
pendix IS] the Sf gyrokinetic equation becomes 



dt 



-a 

"''^ dw\\ 



ZsefMsd{x)s f ^ 

jMsV^ 



T, 



dt 



Ts dip 



Ts dip 



1 dus mslw\\ dfli^ 
Vslhp ^ BTs Ihp 

nisw'^ 3\ 1 dT. 



2T, 2j Ts dip 



.(19) 



Here, the turbulent components of the electromagnetic 
field are contained in 



The ip and a components of the drift velocity are 



'-'d,s 



cl d^o ^(^N + dB 



B de 



BQ, 



2BmQW\\dR IR^ldR 



de 



c : 

iis de 



de 



b- V6I 



and 



47rw;^ dp 

" BQs di) 
mscRQ'^ 
Z.e 



a$o 9$o (Va xh)-Ve 

Ihp ^ ~de B 

dB dB (Va x b) • 

d^ " ~de 



B 



(yR X C) • Va 



dR dR (Va x b) • V6i 

Ihp ~ 'de B 



The parallel acceleration is 

dB Z,,e d<i>o 



de nis do 



RQi 



,dR 

'de 



(25) 



(26) 



b • V0. (27) 



{x)s = Jo{zs) (^<P~ ^w'll^ll^ + 



The Fourier transform of the gyroaveraged collision op- 
erator for collisions between species s and s', {C'jfJ,)s, is 
given in Appendix [C] 



2Ji{zs) nis^ The local versions of the quasineutrality equation ( 14 1 , 

^11' (^^-^ the parallel current equation (15) and the perpendicular 



where J„(zs) is the n-ih Bessel function of the first kind. 
The function Zs{k^, ka, e, /j, t) is 



Zs{kti), ka, e, /i, t) 



kj_^/%iB 

n'. ' 



(21) 



current equation ( 16 1 are obtained in Appendix |Dj and 
are given by 



'^2nZsB j dw\\ dnhsJo{zs) - 



Ts 



0, (28) 



with 

k^ = (fc^lVV'l' + 2k^kaVip ■ Va + A:^ |Vani/2 (22) 

the magnitude of the perpendicular wavenumber. Note 
that the perpendicular wavenumber increases with time 
because Va = . . . — Vip{dft(^/dip)t. This is equivalent to 
the procedure generally employed in local Sf simulations 
to account for the shear in the Ex B drift |Tni[IIJ[57]. The 
generalized potential (x)s appears in the radial turbulent 
drift 



fci^ll ^ ^ ^!L^i^ I dw\\d^hsW\\jQ{zs) (29) 



and 



BB^ 
47r 



271171 sB^ J dw^^d/j hsfJ, 



2Ji{zs) 



0. (30) 



The perpendicular current equation ( 30 1 has been rewrit 



ten as a perpendicular pressure balance. 



= ikac{x)s 



(23) 



and in the nonlinear term 



{{x)s,hs}^c {k[pka - k^k'J{x}s{k'^,k'a) 

xhsik^-k'^,ka-k'J. (24) 



D. Radial flux of toroidal angular momentum 



by 



The radial flux of toroidal angular momentum is given 



(ns)At + (nB)At, 



(31) 
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where 



U,^m,{{R I d^wfsiv-C){^-Vi,)) ) (32) 



/ 4,, 



is the momentum flux due to the species s, and 



1 

An 



i?[(B + (5B) •C][(B + <5B) • V^] 



(33) 

is the momentum flux due to the Maxweh stress. Here, 
(...)^ = {V}-'^ f dedCJ{...) is the flux surface av- 
erage, V = dV/dtp = 2nfd9J, and (. . .)a,a = 
^V'" Va^ '^^ (■ ■ •) and (...)At = At-^ J^^dt{...) are 
coarse grain averages over a radial region Aip and over 
a time interval At larger than the decorrelation length 
and time of the turbulence, respectively. In the local 
approximation, Ai/j is much shorter than the charac- 
teristic length of variation of the background profiles 
of density, temperature and velocity, and At is much 
shorter than the transport time scale. We remind the 
reader that the velocity v = w -I- used in equa- 

tion ( 32 ) and in some other equations in this section, 
is the velocity in the laboratory frame. Note also that 



the distribution function fs includes the lowest order 
Maxwellian /ms, the turbulent contribution Sfs and neo- 
classical higher order pieces. Finally, it is worth pointing 
out that our definition of radial flux of toroidal angular 
momentum does not have units of pressure x length, but 
of pressure x (length)^ x (magnetic field) because our ra- 
dial variable ip, which appears in the gradient V?/", has 
dimensions of (length)^ x (magnetic field). 

To calculate an expression for the turbulent contribu- 
tion to lis, it is convenient to use the full Fokker-Planck 
equation 



^+v.V/. + L,vxb+^ 
ot m„ 



-SE + -w X (5B 
c 



V($_i -^$o) 
V„/, = ^ (34) 



The turbulent Lorcntz force is expressed in terms of the 
peculiar velocity w = v — Rn,(^<^ because the fluctuating 
magnetic fields are calculated in the rotating frame. By 
taking moments of ( 34 1 we obtain a simplified expression 
for [nillS]. The msi?^(v • 0^ moment gives 



2Z.,eR 



m,R^ / d^w /,v(v • CY 



where we have used RB x ^ = Vi/i and V(i?C) = (^-R)C ^ C(V-R)- Taking a flux surface average and the coarse grain 
averages over radius and time in this equation, we find that the second and third terms in the right side vanish to 
order e|psi?|VV'|- The last term only gives a neoclassical contribution, leaving the turbulent transport of momentum 




msR-' I d^w hs{w ■ C + R^c) [SE+ -vfx SB] -C 




Alp/ At 



',PsR\Vi;\. (36) 



Note that only the fiuctuating piece of the distribution 
function 5fs — hg — {ZsE(j)lTs)fMs contributes to the 
turbulent transport of toroidal angular momentum be- 
cause the coarse grain averages satisfy {{S fs) A-tp) At = 0, 
{{SE) A,p) At = and {{SB)A^)At = 0. In addition, 
the contribution from the Maxwell-Boltzmann response 
— {Zse(j)/Ts)fMs vanishes upon integrating over velocity 
space and over the spatial scales of the turbulence. 
The turbulent contribution to the Maxwell stress is 



TTtb 



1 

47r 



R[SB-C)[5B-Vi>) 



(37) 



Ai> 



This contribution is usually small for low /? plasmas, but 
it can become important where the pressure gradient ap- 



proaches the ballooning threshold [29l |30] . 

To study the different contributions to the turbulent 
transport of momentum in ( 36 1 and ( 37 1 , we rewrite 
and using SE = -V_l0, SB = i?||b + V A\\ x b, 

_B||(w X b) and 



SB = 

w X 5B = V_lA||W|| — (wj_ • V_L74||)b 
RC = (Ih/B) - (b X VtP/B), giving 

s{R'[ 



Tjtb „ 
Il„ = TT 



s,|| + T^s,± 



and 



n 



tb 



B = ^B,\\ + TTsa- 



(38) 



(39) 



The transport of momentum has three contributions: 
a piece '^^'ms{Rl^)^^l(^gs^ where gs is very similar to 
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but not exactly the radial flux of particles, a piece 
Ss('''s,ll that is related to the radial-toroidal com- 

ponent of the viscosity and the Reynolds stress, and the 
Maxwell stress tt^h -t-Tr^^j,. The different contributions 
■""sJIj '^s,±, ""-B-II and 7rs,_L are classified according to the 
component of the momentum that they are transporting. 
The components 



and 



-(w X h) ■ 'Vip]{w± ■ V^ii^ 



(:^^ii(v^iixb)-vv. 



Alp 



(40) 



(41) 



A'0 



transport parallel momentum and dominate in a tokamak 
where the magnetic field is mainly toroidal. Here the 
radial component of the turbulent velocity is 



The components 



JivAii xb- Jiw^ ) -W. 



and 



B 



(42) 
(43) 



Alp 



7rB,± = -U^(^h ■ X b) • 

\ ^ ' ''i' I Alp 

(44) 

transport perpendicular momentum and dominate for 
tokamaks where the magnetic field is mainly poloidal. 
Finally, is 



1 



5s = 



(45) 



In AppendixlElwe evaluate tTs,\\i t^s.x^ 9s, t^b,]] and 7rB,± 
employing the iburier analyzed distribution function and 
electromagnetic field obtained from the local formulation 



in section II C In what follows, the different contribu- 



tions to TTsji , TTs.j^ and Qs are classified according to which 
type of fluctuation in the electromagnetic fields drives 
them. The flux of parallel momentum due to species s is 



vr n + TT II 
s.\\ s,\\ 



TT I! , with 
sir 



ka J dO J(j){k^,ka) 



X / dw\\d^hs{-k^,-ka)w\\Joizs), (46) 



V 



X I dw\\dfihs{-k^,~ka) 
and 



^ fca / d6 JA\\{k^,ko 



w^^Joizs) - fJ.B- 



(47) 



Bu Air'^imlcl 

TT = 



1 f 

— ^ fca / d6JB\\{k^,ka) 



Z,eV 



X / dw\\d^ihs(~k^,~ka)'w\\ii — J_L^^ (43) 



Similarly tTs^x^ = 7r^j_ + 7r^\^ + 7rf\^, with 



tt: I = - 



47r^rn^c^ 
ZseV^ 



7- E dejk^4>{k^,k^) 



X I dwiid^j,hs{-k^,-ka)fJ,'^^^^^-^, (49) 



j = E / d9Jk'^Mk^,k^) 

X / dw\\d^hs{—k^,~'ka)'W\\^ — (^59) 

./ Zs 



and 



r k^ j dejk^B\\{k^,k^) 



Z2e2v 



X J dwiidfihs{~k^,-ka)fJ.'^G{zs), (51) 

where G{zs) = [8Ji{zs) + 4:ZsJ2{zs) ~ 4:ZsJq{zs)]/ z^ and 
k"^ — k_L • Vi/' = /c^|V?Ap -t-fccVi/"- Va is the projection of 
the wavevector on the radial direction. The component 

A B 

9s ^9t +9s " +gs ", where 



3s = \nim\ E ^" / d0JR^B(p{k^,ka) 
^ k^, k •' 



X dw\\dfihs{~k^,~ka)Jo{zs), (52) 



47r^z 



E ^« / dejR^BA\\{k^,k„) 



X / d^/is(-fc^, -fca)?!)!! Jo(2;s) (53) 



and 

_B|l 47r^imsC 



E fc^y dejR^BB\\{k^,k^) 



X / dw\\diJLhs{—k^T—ka)lJi — ^-^^.(54) 
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(a) 



(b) 





FIG. 1: Equilibrium (a) and its up-down reflection (b). The 
dash-dot line is the axis of symmetry. 



TABLE I: Equilibrium quantities under up-down reflection 
Before up-down reflection After up-down reflection 



J, Us, $0, P, dp/dtp\R 
B-V9 



B 



J, ris 
B 



,, $0, P, dp/dtlj\R 



TABLE II: Transformation I and II that leave the local 5f 
gyrokinetic system of equations unchanged. 



Before transformation 


After transformation 


Geometry 


Up-down reflected geometry 


I ^5 , 


fl, ka , A\\ 


dQ(/di), k^, 


— n^, —dQ,(/di), — W||, —k^, 


hs, (t>, -B|| 


-hs, -(f), --B|| 


Geometry 


Up-down reflected geometry 


/i, k^, ka, hs, (p, -B| 


/I, k^, ka, hs, 4>, B\\ 




-J, -^1 



The Maxwell stress contribution to the transport of par- 
allel momentum is 

■^B,\\^-^ X! j de JB\\{k^,ka)A\\{-k^,-ka). 

(55) 

Finally, the Maxwell stress contribution to the transport 
of perpendicular momentum is 

^B,± = ^ f dOJk'^\A^^{k^,k^)\^ (56) 



information is contained in Tables [T] and |III| Table III 



III. SYMMETRY OF THE LOCAL 5f 
GYROKINETIC FORMULATION 

In this section we present two symmetries of the local 
Sf gyrokinetic formulation. Both of them involve an up- 
down reflection as represented in Fig. [l] Note that the 
shape of the flux surface is an up-down reflection of the 
original one, but the direction of the poloidal magnetic 
field is unchanged. Thus, this reflection leaves the di- 
rection of the toroidal plasma current unchanged. This 
up-down reflection does not affect the toroidal magnetic 
fleld. In fact, one of the symmetries that we present 
here will leave it unchanged whereas the other reverses it. 
Importantly, for a given equilibrium, with some chosen 
poloidal angle definition 0, it is possible to find a related 
poloidal angle for the up-down reflection such that the 
functions of ip and 9 given in Table |l] satisfy the relations 
given in that Table. We prove this in Appendix [F} 

There are two transformations, which we will call 
transformations I and II, that leave the local 5f gyroki- 
netic system of equations unchanged. They are given 
in Table |ll] To prove that these symmetries hold, it is 
necessary to know how certain coefficients in the gyroki- 
netic equations transform under up-down reflection. This 



derived in Appendix [Gj is for the coefficients that con- 
tain Va. Since they depend on / and dfli^/dijj, they 
change differently under transformation I or II. With the 
information in Tables [l] and III it is tedious but straight- 
forward to show that the formulation given in subsec- 
tion III CI is invariant under the transformations in Ta- 
ble [n] The symmetry under transformation I is the one 
that is closer to the the symmetry discovered in Ref. [15j . 
Since in Ref. [15] some of the effects that we are consid- 
ering here were not treated, the transformation is not 
exactly the same. In particular, Ref. neglected v4|| 
and -B|| and proposed a transformation composed only 
of the up-down reflection and the change w\\ — ?> —w\\. 
This symmetry is only valid for the linearized equations 
and for the particular case = 0, which is the most 
commonly studied in local Sf gyrokinetic linear theory. 
In addition, Ref. 15 did not state that hs and (j) should 
change sign under the transformation. The reason is that 
the symmetry of the linearized equations cannot deter- 
mine if hs, (f), A\\ or B\\ must change sign. Only the 
nonlinear term can determine that. 



IV. CONSEQUENCES FOR THE TRANSPORT 
OF TOROIDAL ANGULAR MOMENTUM 



According to transformation 
hs{k^,ka,9,w\\,^i,t), 
and B\\{kTp,ka,9,t) 
tions in subsection 
then —hs{^k^,ka, 

Aii{-k^,ka,9,t) and (-/c^, /cq, 6*, i) 



I in Table |TT1 if 
ka,9,t), A\\{k^,ka,9,t) 
are solutions of the equa- 
|II C| in a given equilibrium, 

,-wii,n,t), -(t){-k^,ka,9,t), 
are solu- 



tions in the up-down reflection with 17,^ — >■ — and 
d^l^/dip — > —dO^c/dtp- This is important for transport 
of momentum because the radial transport of toroidal 



angular momentum, given in subsection II D changes 
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TABLE III: Changes under transformations I and II. 



Before transformation 



After transformation 



jVaj", (Vq X B; 
I {VR X C) • Va 

Vip ■ Va, V6' ■ Va 



ve, 



|Va|", (Va X B) • V6», 
(V7? X C) • Va 

-Vip ■ Va, -V6' ■ Va 



Vi/) • Va, Vf • Va, |Va|^ Vi/) ■ Va, VS ■ Va, |Val^ 
II (Va X B) ■ Ve (Va X B) • Ve 



(Vi? X C) • Va 



-(VT? X C) • Va 



under transformation I. The contributions tt 



sign 



tj-' ■^f!'!' ^B,\\ and 7r_B,_L change sign. 

The contributions g^, Qs " and " remain the same, but 
since they appear in the momentum flux multiphed by 

their total contribution also changes sign. Thus, for 
any initial condition hg{k^,ka,9,w\\^ fi,t = 0), there is 
another initial condition —hs{—k^,ka,9,—w\\,fi,t = 0) 
that gives the opposite tt., ||, tt^^j^, t^b.\\ and 7rB,i_ and 
the same gs for the up-down reflection with — > — 
and dVl(^/dij} — —dQ.^/dip. Since t^s,\\j ^s,_l, 3s, ''^b,\\ 
and 7rs,± are time averaged over times longer than the 
decorrelation time, they are independent of the initial 
condition, and the radial transport of momentum must 
change sign under up-down reflection and the change 

—7- and dVlc^/dij] —dilQ/dtl'. Because the sym- 
metry changes the sign of each contribution tt^ 

^^I'j' '^t,±' ^fX' ^s.±^ '^B,\\ and 7rB^_L independently, this 
property holds for each one of them independently. The 
same argument is valid for gf, ijs " and " . In the 
particular case of up-down symmetric tokamaks, this 
symmetry implies that the sign of the radial transport 
of momentum changes when the sign of fi^ and (9J7^ /d^p 
is changed, and when fi^ = and dVlc^/dip = 0, it must 
vanish. A consequence of this result is that reversing 
the momentum input in up-down symmetric tokamaks 
exactly reverses the rotation profile Q,q{iP). 

For up-down asymmetric configurations, there is tur- 



TT 1 1 , TT 1 1 , TT 



bulent transport of momentum even for Vt 



C 



and 



dVl(^ I dip = [S] . This momentum flux is strongly depen- 
dent on the temperature and density gradients that drive 
the turbulence. An explicit expression for the turbulent 
momentum flux due to up-down asymmetry has not been 
obtained, but it is possible to comment on the size of this 
contribution for small up-down asymmetry. We distin- 
guish two types of equilibrium quantities that enter in 
the coefficients of the gyrokinetic system of equations in 
subsection II C the quantities Q{ip,9) for which, in an 



up-down symmetric configuration, there exists a poloidal 
angle definition 9 such that Q[tp,9) = Q{'ip,—9), and 
the quantities Q{ip,9), for which Q{'ip,9) = —Q{tjj,—9). 
The former quantities are Q(V', 6*) = {| VtAIMV^P, B, B • 
V6', ris, <I>o,p, and the latter quantities are 

Q{'4>,9) = {Vt/i • V0}. Any equilibrium quantity can 
be divided in its up-down symmetric and antisymmet- 



ric pieces, given by Q^{ip,9) = [Q{ip,9) + Q{iIj,-9)]/2 
and Q^(V',6') = [Qii',9) - Q{ij,~9)]/2 for Q, and by 

Q^{ip,9) ^ [Q{ijj,9) + Q^,-9)]/2 and = 
[Q{tp, 9) — Q{4>, -'9)]/2 for Q. In a tokamak close to up- 
down symmetry, there is a poloidal angle definition 9 for 

which ^ and Q <C Q . In the gyrokinetic equa- 
tion (19) and the Maxwell equations ([28|-(30), the quan- 



tities Q and Q enter in different coefficients. For small 
up-down symmetry, we can expand in /Q^ ^ 1 and 
Q /Q ^1, leaving to lowest order a system of gyroki- 
netic equations that only depend on and Q . In this 
case, the lowest order solutions (p^°\ A''^^^ and b'^^^ 
have the symmetry properties of the up-down symmetric 
case and the momentum flux is zero. For example, 



0(0) _ 4:n^imscl ^-^ 
II — —. / ^ fco 



d9J^(P'^'>\k^,k^) 



X I dw\\dnhf\-k^,-ka)w\\Jo{z^g)=Q, (57) 



where by using the superscript S we have emphasized 
that to this order and Zg are evaluated using only the 
up-down symmetric coefficients and Q . The next 
order correction to the system of gyrokinetic equations 
is linear in the next order corrections h')^\ 4'^^\ A^^^ and 



b\^^ to the lowest order solutions, and contains the anti- 
symmetric coefficients and . As a result, 



.(1) 



A 



(1) 



hf^ 



(1) 

(0) 



< 1. (58) 



Since the lowest order momentum flux is zero, the contri- 
bution due to these higher order correction becomes the 
only contribution. For example, to next order we find 



^' k k 
X j dw\\d^h^g'\—k^ 



V 



J2 / d9J^(P^^\k^,k„) 

ka)wilJo{Zg) 

J2 I d9J^(j)'^''\k^,k^) 



X j dw\\dfih[^\—k^,,—ka)w\\jQ{z^) 
'^^^"^ E Kj d9jU('\k^,k^) 
X I dw\\d^h''g\—k^^~ka)w\\Jo{Zg) 
f d9J^^''\k^,k^) 



V 



X dwiidfMhf\-k^,-ka)wiiJi{z^)zf, (59) 
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where the superscript A in and z-f indicates that (P(k^, kc,, W|| , 0)) 



At 



these are the next order corrections in Q /Q <^ 1 and 
] IQ <C 1 to the coefficients J and Zg, and as a con- 



sequence, they contain Q and Q . Note that tt^;, ' con- 



,0(1) 



tains terms that have exactly the opposite symmetry to 
the symmetry that made the lowest order contribution 
TT^'jl'^ zero. Thus, the contribution to momentum flux 
due to up-down asymmetry is proportional to the asym- 
metric component of the equilibrium coefficients, that is, 

proportional to Q"^ and Q . 

To demonstrate that the symmetry based on transfor- 
mation I is observed in local Sf gyrokinetic simulations, 
the local Sf gyrokinetic code GS2 [3 was used to solve the 
gyrokinetic system of equations ( 19 1 and ([28]) for an up- 
down symmetric equilibrium with $7^ — 0, but with three 
different values of dil^/dip: zero and ■izQ.2q^vti/ R^Bi^. 
Here vu = yjxjmi, Ti and rrii are the ion thermal speed, 
temperature and mass, i?o and Bq are the major radius 
and the magnetic field at the magnetic axis, q = 1.4 is 
the safety factor, and e = 0.18 is the inverse aspect ra- 
tio. All other equilibrium parameters correspond to the 
Cyclone base case [3T]. Only electrostatic fluctuations 
are considered and a modified Boltzmann response |32] 
is employed for the electrons. From the solutions for /i^, 
we calculate the dimensionless function P(k^, ka, 0, W||), 
defined as 

27rcfca0(fc^,fca) f , , f , I, ^ 

P=^^ ,|, / dfih,{-k^,-kc,) 

efRoHiVu (IVV-I)^ J 

X (^^Iw^^Mz,) - ^k^^'^Jl^^ ^ (60) 

where rii, pi — riiTi and pi = miCVu/eB are ion density, 
pressure and thermal gyroradius, and = pi/Ro- The 
function P can be understood as the dimensionless contri- 
bution to the radial momentum flux due to fluctuations 
at the poloidal angle 6 with dimensionless wavelength 
= kjjjei{eR^BQ I q) and k^ = kaei{qle) when inter- 
acting with particles with dimensionless parallel velocity 
W|| = ww/vti- The function P is related to the turbulent 
toroidal angular momentum flux by 



II 



^t_L)At 



27r 
V 



dOJ 



dw||P. (61) 



From the symmetry under transformation I in 
Table |ll] we flnd that (P(k^, kc W|| , 0))At — 
— (P(— k^, kc, — W||, —9)) At for the case with dil(^/dtp — 0. 
The change in sign in the angle is due to the up-down 
reflection. As proven in Appendix |Fj for every angle 
9 defined for an equilibrium, there is an angle 9 for 
the up-down reflection such that Table |T] is satisfied. 
In this case, for which the up-down reflection is the 
flux surface itself, we choose the poloidal angle def- 
inition 9 so that the corresponding poloidal angle 
for the up-down reflection is simply —9. The result 



-(P(-k^,kQ,,-W||,-6'))At is 
confirmed by Figs. [2]ja) and [2]jb), which clearly show 
that the integrals of the function P(k^, k^, 6*, W||) 



Pw||,e(w||,( 



and 



2n 
V' 



d0J / rfwiiP 



At 



satisfy the relations 

Pw||,e(w||,e') = -Pw||,0(-W||, 

and 



k,/; , k, 



(k0,k„) = -Py ^ 



(62) 



(63) 



(64) 



(65) 



This symmetry is broken for the cases with fi- 
nite d^^/dip^ which are shown in Figs. [2]jc)-[2|f ) . 
Observe that the symmetry also gives that the 
case with negative dVl^/d^ has P- (k^, k^, W|| , 6*) = 
-P+(-k^, kQ,-W||,-6'), where P+(k^, kc,, W|| , 6*) is the 
function P for the case with positive dil(^/di/j, and 



P_(k^, ka, W||, 0) is the function P for the negative 
d^Q/dijj. This is confirmed by comparing Fig. [2][c) 
with [2]je), and Fig. [2]jd) with[2]jf). As a consequence 
of these symmetries, the case with dilQ/dip = has 
(tt^II -I- 7rf_|_)Af = 0, while the other cases have {nf^^ + 

<i)At = ±2.6efi^oP^(|V^|)^. 

Following the same reasoning as for transformation I, 
transformation II in Table [ll| implies that the transport of 
momentum remains the same for an up-down reflection 
and a reversal of the toroidal magnetic fleld B/^ = I/R. 
This means then that for an up-down symmetric toka- 
mak, the transport of momentum cannot depend on the 
direction of the toroidal magnetic field. For example, the 
Coriolis pinch [7J will remain the same under a reversal 
of the toroidal magnetic field. 

These symmetries of the lowest order gyrokinetic for- 
mulation have very important consequences for the mod- 
eling of spontaneous rotation profiles in tokamaks with- 
out net momentum input. When no momentum is in- 
jected, the radial momentum flux across any flux surface 
within the separatrix must vanish, i.e., = 0. In the 

absence of up-down asymmetry, the radial transport of 
momentum given by the lowest order gyrokinetic model 
vanishes for Jl^; = and d^lc^/di/j = 0. Thus, for zero or 
negligible velocity in the edge, the solution to the equa- 
tion Il{il(^,dfl(^/dip) = is ^({i/j) = everywhere. To 
obtain ri^(?/') 7^ with the lowest order gyrokinetic for- 
mulation, we need to assume that there is some up-down 
asymmetry [5] . This limits the applicability of the lowest 
order formulation to studies of spontaneous rotation in 
very up-down asymmetric configurations. When higher 
order terms are considered |17[ I18| , more mechanisms for 
the generation of spontaneous rotation arise. 
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FIG. 2: (Color online) Contributions to momentum transport for difTerent values of dQ.Q/dip: (a)-(b) d^lt^/dip = 0; (c)-(d) 
dQ.c^/dip — 0.2q^vti/e^RoBo; and (e)-(f) dfl^/dip = — 0.2 q^vu/e^ Rq Bo ■ Note that in figures (e) and (f) the color map has been 
inverted to aid the comparison with figures (c) and (d). 
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Appendix A: Collision operator 



The linearized collision operator is given by 



(Al) 



The different contributions to the collision operator are 
the pitch-angle scattering 



the energy scattering 



(A2) 



777/ TJJ 

TL' = ^ss' w ( w • V^hs + ) , (A3) 



and the integral contribution 



T, 



It 



(A4) 

where the Rosenbluth potentials Hs' and Ls> are obtained 
by solving the differential equations Wf^Ls' ~ —8TThs> 
and VlHs' = Ls'. Here -/ss' = 2nZ^Zle'^\nA/ml, In A 
is Coulomb's logarithm, 



23/22-3/2 I 2^5, ^ ^ 2wt, dWs 



is the pitch-angle scattering collision frequency, and 



23/22^3/2 1 y^5^ ^y^^ 



is the energy scattering collision frequency. We have 
used the abbreviated notation Wgi = w j y/lTgi jra^' ■, 

= (2/^/7^) /J"'' dy exp{-y'^) and d<^ldws' = 



Appendix B: Derivation of the local 5f kinetic 
equation ( |19[ ) 



In this Appendix we derive (19) from m\. We Fourier 



analyze the equation term by term. 

We change from the spatial coordinates X to the new 
coordinates Q;(X,t) and 6'(X). After this co- 

ordinate transformation, and after Fourier transform- 
ing, the term dhg/dt + C ' Vx^s becomes dhs/dt + 
[{da/dt) + RQ,^C, ■ \/a\ikahs. We drop the subscript x 
in Vx for simplicity. Employing da/dt + RVli^C, ■ Vq = 
0, dhs/dt + RVli^C, ■ Vyjis becomes dhg/dt. Similarly, 
the Fourier transform of d{x)s/dt + R^1(^C ■ Vx(x)s is 
d{x)s/dt. 



The Fourier transform of Vd,s- Vx^is is ifc^/isVd,^- V-^-l- 
ikahsVd.s ■ Va -|- Wd.s ■ y9{dhs/d9). The contribution 
Vd,s • V6{d/d9) is negligible since we assume that the 
parallel derivative d/dO is small compared to k±^ . Before 
evaluating v^^^g • Vip and Wd,s • Va, we rewrite the cur- 
vature drift by employing the MHD equilibrium result 
b-Vb = B-^Vi^B+{ATT/B'^)(dp/dii\R)Vii. Then, defin- 
ing = Wd,s ■ V?/' and v"^ ,, = v^^s ■ Va, and employing 
that the gradient of the axisymmetric quantities <i>o, B 
and R can be written as Vx = \7^p{d/d'i(j)+'^0{d/d9), we 
find expressions (251 and ( 26 ) . To obtain the final form of 



these equations we have used (V^ x ^) • V'0 — i?B • V6 for 
the ip component of Vco.s, b x Wijj ~ Ih~ RBC, for the ?/; 
component of v^o, vm,s and Vc/,s, and (VQ;xV?A)-b = B 
for the a component of ^eo, ^m,s and "Vcf.s- 

The nonlinear term {(x)sj ^s} is the Fourier transform 



of s ■ Vx^s- The generalized potential {x)s in (10 1 can 
be Fourier transformed by realizing that 



{X)s = (Xexp(i5(x))), 



(Bl) 



where x — 4>^c "'"(ti;||yl||+wj_-Aj_) and 5'(x) = k^ijj{-x.) + 
fc„a(x). Using S'(x) = S[~K + Ps) ~ S'(X) -f k_L • p., with 
kj^ = fc^VV' + fcjjVa, we find that the Fourier transform 



of (Bll is 



27r 



27r 



) / dif exp(zkj_ ■ Ps 



— Aj^ • / wj^ exp(ikj^ • Ps)- (B2) 

TTC J 



To obtain the final result in (20 1, we need to employ 



B 



II 



and 



ih ■ (k_L X A_l), 



j di^ exp(ik_L • Ps) = 27rJo(zs) 



(B3) 



, ,., \ o 'i-nT'sCfj.2Ji{zs) 
d(fw± cxp(ikj^ • Ps) = — 27r^j b x kj^, 



(B4) 



where Zs is defined in ([2l]). To prove (B3|, we use (11) 
to write k_L • = -{^/2jlB /Qsl'^^diik^'fpl'^ + ka'^a 
VV') sin if + Bka cos </?]. Then, 



k_L • Ps = -Zs s'm{ip + ipk), 



(B5) 



with tpi^ = arctan[i3fcQ,/(fc^|V'(/'p -I- k^Va ■ V?/;)]. Using 
this expression and 

oo 

exp(izsiniy9) = J„(z) exp(m(^), (B6) 



n— — oo 



we obtain (B3). Equation (B4| is obtained by realizing 



that Wj_ exp(zkj^ • ps) — iilsb x Vk^ exp(ik_L • Ps) and 
hence J dLpyv±^ exp(ikj^ • ps) = 27rir2sb x Vk^ Jo(-Zs)- 
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Finally, the last term in the right side of ( 19 ) is the 
Fourier transform of the last term in ([s]). To obtain the 
final result in (19), we use ^ to find 



1 dus Zscd'Pa msRVildR 



30 



+ 



n 86 



86 



= 0. 



(B7) 



Appendix C: Gyroaveraged collision operator 

In this Appendix we Fourier transform the gyroaverage 
of the collision operator given in Appendix]^ To obtain 
the Fourier transform, we need to express the collision 
operator as 



(CI) 



where S'(X) — /c^;-i/;(X) + kaCnilC) depends on the gy- 
rocenter position X. This Appendix contains the pro- 
cedures that we use to write the different terms in the 



gyroaveraged collision operator as in (CI I. Once the op- 
erator is expressed in this way, it is easy to see that the 
coefficients {C^^],)s are the ones that appear in 

To differentiate with respect to w in Vu, • F^, and 
• f , we use that 



/is = ^ /is exp(i5(X)). 



(C2) 



Employing "S/^w^ = b, V^/i = w±/B and Vu,S'(X) — 

VujX • kj^ = il^^b X kj^, we find that the Fourier trans- 
forms of the gyroaveraged divergences of the vectors f 
and Tig, are [55] 

(V. • r^,)s = + - '-^Di, (C3) 



8wii 



8^i 



„ ^ / 8hs „ 8hs ■nis'up , 



' 8w\\ 



8ti T,, 



(C8) 



and Dig, = i^gg,hs. Note that ^ w^^ + 2fiB. 

To Fourier transform the gyroaveraged divergence of 
r^, , it is convenient to write the Rosenbluth poten- 
tials Hs' (x, w\\, fi, if, t) and Ls' (x, , /z, (f, t) as functions 
of the position of the particle x instead of the guid- 
ing center X. To obtain Hg, and Ls' fro m eq uations 
y^Ls' — —SttJis' and S/^Hg' — Ls', we use (C2|, where 



5(X ) — S{x. — ps) ~ 5(x) — kj^ • ps, and wc employ (B5) 
and ( B6 ) to write 



hs=^ ^ Jn{zs)hsexp{in{ip+Lpi^)+iS{x)). (C9) 

,ka n=~oo 

Using the decompositions 

oo 

Ls' = ^ ^ Ls',nexp{in{(p + (fii^) + iS{x)) (CIO) 



fcj, .fc„ n=—oc 



and 



Hs' ^ E Hs',nexp{iniip + ipi,)+iS{x)), (Cll) 

,ka n= — oo 

we obtain the functions Ls'^n{k^,ka,w\\, fi,t) and 
Hs> ,nikif,j kajW^^, ^,t) from equations 



8^Ls',n 2__9_/ 8Ls'.n \ 



and 



-8T:Jn{Zs')hs' 

(C12) 



and 



II 



Of, n^s 



Here the pitch-angle scattering terms are 



dhs dhs 



2^^'^ss'[B^-w 



117577 ' 



^11 



dhs ^ ^\\dhs 



(C5) 



(C6) 



With the decompositions (ClOl and (Cll I, and Vu,W|| 



b, V„Ai = w_l/B, Vy,ip = -w X h/2fj,B, VyjV^w\\ = 0, 

^w'^wfJ- ~ (I -bb)/_B and V^Vu_,<^ = [w_L(vkf x b) -|- 
(w X b)wj_]/4^2^2^ obtain 



v.-e= E E 

fcj, .fee r! = — oo 



ss' ,iL)^ I ,n 



8C^ 



exp(m(c/j + 9?k) + iS'(x)),(C14) 



where 



and D^g, = [1 -I- {w'^^/ fiB)]i'^g,hs, and the energy scatter- 
ing terms are 

/ c)h dh in w'^ \ 

rLv.„ =^ii-L'Ui^ + 2a.^ + ^/is 1, (C7) 



^ ss' I ,n 



rrisjss'fMs / d'^Hs'^n 



Ts 



w 



II" dwf^ 



^ s' .71 . J- s ^^s' ,n 



(C15) 
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djj? 



BT, 



Appendix D: Derivation of the local Sf equations for 
the electromagnetic fields 



and 



z2 ^BT, 



dfj, nis' dfi 



(C16) 



In this Appendix we show how to obtain ( 28 )- ( 30 1 from 



dw\\ 



dHs 



(C17) 



Using S'(x) = ^(X + ps) ~ S'(X) + k_L • and ([B6]), we 
find 



(14|-(16). We only need to Fourier transform (14|-(16|. 
In equation (14 1, the integral / dw\ \ dfid^phg can be writ- 
ten as 

j dw^^dfihg j exp(zS'(X)), (Dl) 

where S'(X) = k^i!{X) + fcc.a(X). Using S'(X) = 
5(x — ps) — S{x) — li± ■ ps and (B3), we find 
J d(p exp{iS{X)) = 27rexp(iS'(x))Jo(z,.), and ^ im- 
mediately follows. Equations (29 1 and (30) are ob- 



taine d f rom ( 15 ) and (161 in a similar way. In equa- 
tion ( 16 ) we find the integral / dip wj^ exp(z5(X)). Using 



exp(iS'(x)) = J„(z^)exp(-m((p-|-^k) + i5'(X)), 5(X) - S{x) - ik± ■ p, and pi]), we obtain ([30l 



n— — OO 



(C18) 

giving that the Fourier transform of the gyroaverage of 



V. • r«, is 



OO / QQR 

n— — OO \ 



(C19) 



This equation combined with (C3| and (C4| gives that 
the Fourier transform of the col 



C. 



lision operator is 
R 



, 1 = — fr^ , + , -I- F", 

/ s dw\\ V '""11 '•'"ii 



Q2 \ ss' ^ ^ss' 



Di 



where 



ss' ,n ~ ^ ^ J'n,{Zs)G ggf 



(C20) 

(C21) 
(C22) 



and 



Lt „„/ 

+ Jn{Zs)Af, 



s' .wi I ,n 



(C23) 



To obtain the final result in (C20), we integrated 
the first two terms in (C19) by parts in w^^ and /i, 



and we used dJn{zs)/d'w\\ = and dJn{zs) / dfi, = 
{kj_^/B /ilsy/2fi){dJn/dzs), where dJn/dz^ = [Jn-iizs)- 
J„+i(z,)]/2. 



Appendix E: Derivation of the local 5f flux of 
toroidal angular momentum 



In this Appendix we show how to obtain ( 46 )- ( 56 1 from 



(40 1- (45). We only show in detail how to obtain (|46|) 



from the term proportional to — (c/B)V(/) x b in (W0| 
We then briefly cover the other terms. Writing the flux 
surface average explicitly, we obtain 

J 



'II ~ V'Aifj 

j^a 



di: dO dC 



B 



xika(l){k,^,ka) J d^whs{k'^,k'^)w\\ 

X exp(if7(x, X, (El) 

where we employ the abbreviated notation 
ct(x,X, fc,0, fca, fc^, fc^) k,^ijj{x) + kaa{x) + k'^'ip(X) + 
fc|jQ;(X). Using X ~ X — Ps, a can be simplified to 
CT ~ (A:^ -I- k'^)il;{x) + (fc^ + fc^)a(x) - ik^ • ps, with 
= k'^Vip + k'^Va. Using a = - qd - n^t, and 
noting that the integral over ip and C is over distances 
much longer than the perpendicular correlation length, 
we find that only the contributions with k'^ = —k^ and 
k'^ — —ka do not vanish, leading to 



2TiimsCl 
V' 



E y d6 j^(l){k^,ka) 



J d^whs{~k^, ~ka)w\\ exp(ik_L • Ps)- (E2) 



To obtain the final result in ( 46 1 , we use d^ 



B dw\ I dfJL dip and ( B3 ) . 

Equations (47)-(SSl are obtained in a similar way. To 



obtain the final expressions, we only need (B3), (B4| and 



2TTflB 



(i(y5 w_l(w X b)exp(zk_L • Ps 



^bx I - — G(zs)(b X k_L)k_L 



,(E3) 
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where G{zs) = [8Ji{zs)+AzsJ2{zs)-4:ZsJo{Zs)]/z^,. This 
result is obtained by realizing that w_l(w x b) exp(ikj^ • 
Ps) = -f^sbx Vk_LVk_L exp(ik_L-Ps), giving / d(^w_L(wx 
b)exp(ik_L • Ps) = -27rrj2b x Vk^ Vk^ Jo(2;s). 

Appendix F: Up-down reflection of the equilibrium 



Appendix G: Proof of the results in Table [ITT] 



In this Appendix we evaluate the coefficients in Ta- 



ble III and express them as functions of other coefficients 
whose behavior under the transformations in Table ITTl we 
know. 

Using Va X B = {Vtp ■ Va)Va - |VapVV', we find 



In this Appendix we show that it is possible to find a 
poloidal angle definition 9 for the up-down refiection of 
an equilibrium such that Table |l] is satisfied. 

We first calculate the geometrical coefficients B, B- V^, 
J, IVV'P, • V6' and jV^p as functions of ip and for 
some choice of poloidal angle definition 9. We invert 
the functions iIj{R,Z) and 9{R,Z), where R and Z arc 
the typical radial and axial coordinates satisfying Vi? x 
S/Z = C. We then write all the coefficients using the 
functions R{ip,9) and Z{il;,9). The poloidal component 
of the magnetic field B • V9 = (VC x V^) ■ V9 is written 
in terms of R and Z as 



B-V9 



R 



dR dZ dZ dR 
d^~d9 ~ IhjjlM) 



(Fl) 



The Jacobian J ^ \B ■ V6i|-i. The products jV^Z-p, 
• V9 and |V6'p are 



dR 
1)9 



dZ\ 



OR dR dZ dZ 



and 



dR 



(F2) 



(F3) 



(F4) 



{P 



The magnitude of the magnetic field is 

For an equilibrium given by R{%p,9) and Z{ip,9), the 
equilibrium described by R{'ip^ 9) and —Z{ip, 9) is the up- 
down reflection. This is a way of choosing the poloidal 
angle definition 9 for the up-down reflection. With this 
choice, equivalent to replacing Z by —Z in (F1)-(F4), 
we obtain most of the results in Table |l] The poloidally 
varying quantities rij, p and dp/dip\ji remain the same 
because they depend on 9 only through the function 
R{ijj,9), which is unchanged upon up-down refiection. 



(Va X B) • (Va • VV')(Va • V9) - \Va\^{Vip ■ V9). 

(Gl) 

Using Vi? X C = -VZ, we find 



BZ dZ 
(VRx C) ■Va = - — Wt/j -Wa- —W9 -Va. (G2) 

dtp o9 



Finally, the coefficients |Vap, Va • V^" and Va • V0 are 



IVaH = 



i?2 



21 



?9- 



t- 



i?2B -^79 \ dip ^dtp dip 

I 



i?2B • V6' 



\/ip-\/9 



\V9\\ (G3) 



Va-Vip 



dq di9 

dip dip ' dip 



i?2B . \/9 



Vijj ■ S/9 (G4) 



and 



Va • V9 



^dq_ (M ^ dn^ 
dip dip dijj 
I 



i?2B • V9 



Wip ■ V9 



\V9\\ (G5) 



with the function d{ip,9) defined in (18) 



With these expressions, and employing Table|T] we can 
show that transformation I, which is an up-down reflec- 
tion and a sign change in fJ^ and dil,(^/dip, transforms 
the coefficients as in Table |III| We can prove the same 
for transformation II, which is an up-down reflection and 
a sign change in /. 
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